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Abstract 

The structure of globular clusters and elliptical galaxies are described in an unified way through 
a new class of lowered models inspired on the nonextensive kinetic theory. These power law models 
are specified by a single parameter q which quantifies to what extent they depart from the class of 
lowered stellar distributions discussed by Michie and King. For ^ = 1 the Michie-King profiles are 
recovered. However, for q < 1 there is a gradual modification in the shape of the density profiles 
which depends on the degree of tidal damage imposed on the model, thereby also providing a good 
fit for globular clusters. It is also shown that a subclass of these models, those with a deeper 
potential and q slightly less than unity, present a distribution resembling the de Vaucoulers r^/^ 
profile which yields a good description of the structure of elliptical galaxies. This subset of models 
follows this trend, with a slight departure over nearly 10 orders of magnitudes. 
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I. INTRODUCTION 



It is widely known that the Maxwell-Boltzmann isothermal distribution is somewhat 
truncated by the galactic tidal field responsible for the capture of stars that become loosely 
bound to the main structure. This means that the ideal distribution function (DF) is 

n 

depopulated due to the stars escaping from the system. As shown by Michie to the case of 



non-rotating, spherically symmetric systems, the resulting distribution may be obtained 
solving the Fokker-Planck equation. This class of distributions were rediscussed by King 
who showed that they give a good representation for the observed density profi 
clusters (GC). The phase space density for such models may be expressed as 
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f{e) = Po(27ra2)-3/2 [e- - l] e<0 (1) 

f{e) =0 £ > 

where the dimensionless quantity e = + $(r))/cr^ is the energy of any given star in 
units of the average thermal energy ($ is the gravitational potential and a is the velocity 
dispersion). Actually, in the case of GC, some independent numerical analysis have 
confirmed the quantitative support provided by such distributions usually called lowered 
models in the hterature jsl. 

Unfortunately, these Michie-King models do not describe many elliptical galaxies whose 
smoothness and symmetry are even more striking than those presented by GC. Actually, 
instead of a specific lowered distribution, the profiles of some giant elliptical galaxies are 
very well described by the de Vaucouleurs r^/^ law which provides sometimes a remarkable 
fit over 10 orders of magnitudes of surface brightness (see, for instance, de Vaucouleurs et 
al. [q^ and references therein). Conversely, there are also some elliptical galaxies (like NGC 
4472) that are not fitted by the de Vaucouleurs law as they are by the lowered Michie- 
King models. Such a state of affairs is quite uncomfortable both from a methodological 
and physical viewpoints. In particular, it opens a large window for adjusting mechanisms 
whose primary objective is somewhat explain the existence of individual profiles even for 
the spherically symmetrical case. 

On the other hand, an increasing attention has been paid to possible nonex- 



tensive effects in the fields of kinetic theory and statistical mechanics 10|(see 



http://tsallis.cat.cbpf.br/biblio.htm for a regularly updated bibliography on this subject). 
The main motivation is the lack of a comprehensive treatment including gravitational and 
Coulombian fields, or more generally, any long range interaction for which the assumed 
additivity of the entropy present in the standard approach is not valid. Inspired on such 
problems, Tsallis proposed the following g-parameterized nonextensive entropic expression 



^. = fei^ ^\^^;f^ , (2) 

where A;^ is the standard Boltzmann constant, Pi is the probability of the i-th microstate, 
and g is a parameter quantifying the degree of nonextensivity. This expression has been 
introduced in order to extend the applicability of statistical mechanics to system with long 
range interactions and has the standard Gibbs-Jaynes-Shannon entropy as a particular lim- 
iting case (q = 1), namely 

S = -ks^Pilnpi . (3) 
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Ten years ago, the first attempts for exploring the kinetic route associated to Tsallis 
entropy approach appeared in the literature. The original kinetic derivation advanced by 
Maxwell (12j] was generalized to include power law distributions as required by this enlarged 
framework isl (see also jl^). In particular, it was shown that the equilibrium velocity 
g-distribution is given by 



1 - (1 - g) 



1/(1-9) 



(4) 
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where the Bg is a g-dependent normalization constant which reduces to the Maxwellian 
value for g = 1. As shown in [l^, the above distribution is uniquely determined from 
two simple requirements: (i) isotropy of the velocity space, and (ii) a suitable nonextensive 
generalization of the Maxwell factorizability condition, or equivalently, the assumption that 
-^('^) 7^ f{vx)f{vy)f{vz). More recently, the kinetic foundations of the above distribution 
were investigated in a deeper level through the generalized Boltzmann's transport equation 
which incorporated the nonextensive effects using two different ingredients. First, a new 
functional form to the kinetic local gas entropy, and, second, a nonfactorizable distribu- 
tion function for the colliding pairs of particles whose physical meaning is quite clear: the 
Boltzmann chaos molecular hypothesis is not valid in this extended framework. It was also 



shown that the kinetic version of the Tsalhs entropy satisfies an ifg-theorem 15j, and still 
more important, the g-parameterized class of power law velocity distributions emerged as 
the unique nonextensive solution describing the equilibrium states. Several physical conse- 
quences (in different branches) of the nonextensive kinetic theory have been investigated in 
the literature, which includes its influence on the transport coefficients (spatially inhomoge- 
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neous dilute gas) of real gases and plasmas [IJ, Il6[ . In particular, a reasonable indication 
for a non-Maxwellian velocity power-law distribution from plasma experiments have been 
demonstrated by Liu et al. 17| and by Lima et al. |l8l| . In high energy physics, many 



studies are illustrated by the solar neutrino problem [ISj], scattering processes in e^e~ an- 
nihilation {20!, heavy ion collisions 21 1, and the charm quark dynamics for a thermalized 
quark-gluon plasma 22], among others (see lo| for an extensive list of applications). Even 
experimental and theoretical deductions of the nonextensive q-parameter have been dis- 
cussed in the literature by taking into account possible dynamical correlations or even the 
presence of inhomogeneiti es ( due to clustering phenomena) thereby driving the system to a 



kind preequilibrium state 24|]. In this sense, the q-parameter cannot be thought as a fltting 
parameter, but a known quantity that should be ultimately determined by the correlations. 
The main lesson is that for systems with sufficiently complex dynamics (for instance, due to 
long range forces) other than Boltzmann-Gibbs statistics can provide a better description. 
In the astrophysical context where our interest is focused, Tsallis' statistics has been 



applied to stellar coUisionless systems [25|, l26| , peculiar velocity 



|27l | , in studies of Jeans instability 



'unction of galaxies clusters 



29j and the Chandrasekhar 



gravothermal instability 
condition for the equilibrium and stability of a star 30|. There are also some expectations 
ihat such an approach may put some light on the universal structure of dark matter halos 
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In a previous paper, we have determined the radial and projected density profiles for 
two large classes of isothermal stellar systems with basis on the equilibrium power law q- 
distributions [3J]. Such models are based in the following phase space density 
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where the g-parameter quantifies the nonadditivity property of the gas entropy. For generic 
values of g 7^ 1, this DF is a power law, whereas for g = 1 it reduces to the standard 
Maxwell-Boltzmann distribution function. Formally, this result follows directly from the 



known identity limd^o(l + dy)d = exp(y) as can be seen in 37]. The quantity Cq is a 
g-dependent normalization constant given by 28 ] 



and 
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which reduce to the expected result in the limit g = 1. 

With basis on the ideas of Michie [l| and King 2|, in this work we propose a new class 
of lowered nonextensive models which are naturally associated with the above equilibrium 
distribution. It will be explicitly assumed that the nonextensive isothermal q-distribution 
discussed in |28| cannot strictly be attained for a real stellar system due to the presence of 
tidal effects with the mean local gravitational field, stellar encounters or any other relaxation 
mechanism. This means that such mechanisms must gradually modify the ideal velocity 
distribution in such a way that the final distribution drop to zero at a finite velocity. Hence, 
one may suppose that the stellar system approaches as far as it can to a quasi power law 
final state which collectively can be termed nonextensive lowered spheres. As we shall see, 
this new class of models may potentially describe both classes of spherically symmetric 
systems, namely globular clusters and elliptical galaxies, thereby reinforcing the possibility 
that nonextensive effects may have a considerable importance in the astrophysical domain. 

The paper is organized as follows. Next section we set up the basic equations defining 
the lowered nonextensive models. In section III we discuss a large class of density profiles 
obtained trough a numerical solution of the Poisson equation. Although reproducing the 
Michie-King models exactly for g — > 1, we show that for g smaller than unity the standard 
profiles are moderately modified, and therefore, they also provide a good representation 
for the structure of GC. In section IV we discuss a specific q-lowered model characterized 
by a relatively deep central potential. It resembles the de Vaucouleurs law being therefore 
potentially useful to describe the structure of elliptical galaxies. This model has an extra 
bonus: it predicts the existence of small fluctuation similar as that ones appearing in the 
profiles of early type galaxies, as discussed by Caon 38] based on the Sersic r^/" law. Finally, 



a summary and the main conclusions are presented in section V. 



II. LOWERED STELLAR DISTRIBUTION: NONEXTENSIVE APPROACH 



Let us now consider a class of lowered nonextensive q-distribution defined by 
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which arises naturally from our earlier adoption of the isothermal power law distribution 
to the unperturbed coUisionless system (cf. equations (1) and (2)). This means that to 
the accuracy of the Michie-King approximation, the above distribution corresponds to the 
simplest non-extensive extension representing the steady state solution of the Fokker-Planck 
equation [35j (see also [36|). The corresponding mass density profile can be obtained by 



integrating this distribution over the whole possible energy interval. It thus follows that the 
density profile is given by the expression 
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where = ^/cr^ is the dimensionless potential, and the integration limits have been defined 
by keeping in mind that only bounded objects are present in the stellar distribution. Once 
again, this expression asymptotically approaches the lowered Michie-King models when q 
1. This expression linking the mass density to the gravitational potential does not have 
a closed analytical form and therefore must be numerically integrated. Now, in order to 
simplify the numerical algorithm, it proves convenient to introduce the variable, e = x0, 
with (Q assuming the form 
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which is more amenable to be solved by numerical discretization. Another benefit from this 
expression is that we can easily verify that the relation between the mass density and the 
gravitational potential has two asymptotic regimes. In the first limit, when ~ 0, we are 



close to the external surface and we can expand the power expression appearing inside the 
brackets using the binomial expansion to obtain 



p^^(-0p 0^0 (11) 

showing that the external structure of these models correspond to a polytropic sphere with 
Lane-Emden index n = 5/2. In fact for —(f) > 1 this expression gives a good representation of 
the model as we can verify from figure [H In the internal structure we can find another regime 
in those regions where the gravitational potential is sufficiently large so that (1 — q)4>x » 1 
in which case, 

poc(-0)^^ (12) 

showing that in this limit the internal structure is also described by a polytropic sphere whose 
Lane-Emden index, n = (5 — 3g)/2(l — g), depends on the nonextensive parameter. In fact as 



showed in 3J], the nonextensive isothermal distributions have a profile density corresponding 
exactly to this polytropic structure. Therefore, in those models were the central potential 
is sufficiently deep, the internal structure is closely described by this approximation. In 
particular, we can verify that the model q = 5/7 corresponds to a polytropic sphere with 
n = 5 which is the limiting case dividing the Lane-Emden family in a branch having finite 
mass and radius from those having infinite radius and mass. Therefore, in the absence of 
tidal truncation, models satisfying the restriction 5/7 < g < 1 have infinite mass and infinite 
radius and are the ones that we will give more attention bellow. We can also see that when 
q = the whole structure is exactly represented by a n = 5/2 polytropic sphere. In figure 
[H we show dependence of the mass density as a function of the gravitational potential for 
a set of representative values of the q index. We can see from this plot that in the general 
case the whole structure can be described as the result of a smooth transition between 
n = 5/2 at the external region to an internal region having n = (5 — 3g)/2(l — q). As 
happens in the Michie-King models, the value of the central potential is also arbitrarily set 
to define implicitly the total mass and the external radius. Therefore, only objects with a 
sufficiently deep gravitational potential have the ability to show the very internal limiting 
polytropic structure. In all the other cases the internal region falls in the transition regime 
where the external polytropic is gently changing towards the asymptotic internal limit. In 



certain sense, this characteristic is responsible by the variety of density profiles found in 
these models. 



III. THE INTEGRATED DENSITY PROFILES 

Having obtained the density as a function of the gravitational potential we can proceed 
by solving the Poisson equation which provides our dynamical link to obtain the density 
profile. Following King [2i] we normalize the density to its central value po and introduce a 
core radius defined as 

9cr2 

(13) 



that will be used as a unit to measure the radial distance. In the normalized radial coordi- 
nate, the Poisson equation becomes 

dr dr 

where the dimensionless density q = p{(j)) / p(0o) can directly be obtained from the solution 
of equation ffTOl) for each choice of the parameter (pQ, whose value will fix the whole structure 
of each model. 

The Poisson equation can numerically be solved by imposing the boundary conditions 
= 00 a^iid d(f)/dr = a.t the origin r = 0. The solution proceed from the center and stop 
at the external surface where r = rt and (f) = g = 0. The external and the core radii define 
naturally a concentration index C = rt/rc indicating how important is the truncation effect 
due to tidal field. Models with low values of C correspond to structures where the tidal field 
had imposed a severe damage to the unperturbed structure. In figure [2] we present a set of 
solutions for three values of q. In the upper panel we show the q = 1 model, corresponding 
to the usual class of lowered Michie-King models. In the other two panels we present the 
profiles for q = 0.9 and q = 0.8, showing that models with the same extent of the standard 
lowered models have different density profiles. Note that all curves have zero gradient for 
sufficiently small values of r regardless of the value of q. Indeed, such a behavior is more 
noticeable for smaller values of q. In figure [3] we present the projected density profile for each 
model presented in figure [21 Assuming that light traces mass these profiles should directly 
be compared with the surface brightness profiles. 
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FIG. 1: The density profiles as a function of the gravitational potential has two regimes. In the 
external region the structure is approximately described by a n = 5/2 polytropic index for all 
models. In the central region of those objects having a deep potential, the limiting structure 
depends explicitly on the nonextensive parameter, and resembles a polytropic with n = (5 — 



The total mass of the q-lowered models can be evaluated by integrating the density profile 

as 



where is a dimensionless mass indicator which can be obtained trough numerical integra- 
tion. Prom Tables I and II we can see that models with lower concentration C representing 
severely damaged structures tend to show a mass distribution close to an uniform spheres 
with /i = 1/3. On the contrary, in objects with higher concentration index the mass is 
heavily dependent on the central core distribution. In the course of the numerical solution 
the integration begins by defining an arbitrary value for the central potential 0o and car- 
rying out the integration to a null value at the external surface. However, since we have 
estimated the total mass we may use its value for determining the gravitational potential at 
the external surface —GM/rt which can be expressed in our units as 



3g)/2(l-g). 
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FIG. 2: Mass density profiles for three representative models. The upper panel correspond to a 
q = 1 model which is exactly the standard Michie-King models. The other two panels correspond 
to models with q = 0.9 and 0.8 respectively. We can see that models having the same size in units 
of Vc, but different values of q, show a large variety of density profiles. In each panel it is also 
indicated the corresponding — (^o of each model. 



(t>t = $(rt)/(T2 = -9///C (16) 

and this value can be used as a correction to be added in order to find the gravitational 
potential at the center. As a final step, the corrected potential can be used to estimate the 
potential energy 

U = - A-K^prHr = (17) 
2 Jo ft 

where 

c 

y = / (f)r'^dr. (18) 



In Tables I and II we show the major parameters obtained in this class of lowered models. 
The first column present the uncorrected central potential, used to start the integration 
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TABLE I: 
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TABLE II: Resume of the lowered model parameters. The upper panel, g = 1, shows the same 
solution as founded with the King models. The other panels illustrate the effect of changing the q 
non extensivity parameter on the structural parameter. 

of the model. The radius containing half the mass, r/i, is shown in column 2 in units of the 
core radius. For the King model, = 1, the half radius is relatively stable even considering 
the large variation covered by the concentration index. We see that the tidal radius, r^, 
varies by more than two orders of magnitudes while the half radius is kept almost fixed at 
70% of the core radius. For the other models this stability of the half radius is not generally 
preserved. As an example, for a model with q = 0.9 a similar variation of the concentration 
index would imply in a variation of a factor of two on the half radius. Actually, models with 
lower q index show a more important departure from the Maxwellian case since the variation 
of the half radius is even larger. This large variation in the half radius is a consequence of a 
change in the shape of the density profile for these models as can be seen from figures [T] and 
[21 In fact, these figures show that for a similar external radius the gradient in the density 



profile is different for eacli value of q. 

In the last two columns of Tables I and 11 we present two concentration indices based 
only on the projected surface density distribution. They are based in the radii containing 
respectively 1/4, 1/2 and 3/4 of the total mass in each model. The concentration index 
Ci is defined to be equal to ri/2/ri/4 and therefore measure the degree of concentration in 
the central mass distribution. On the other hand the index C2 is estimated as r^n/rii2 
and measure the spread of the mass distribution of the outer region. These two indices 
are useful to give a broad representation of the gamma of variations seen in the surface 
density profiles for these models. For a single parameter family each profile is represented 
as one single point in the Ci,C2 plane. This feature is illustrated in figure H] for each model 
represented in Tables 1 and 11. For a given value of q the models falls in a line representing 
the gross features of the individual profiles. Along each line q line the position of the points 
are uniquely determined by the parameter — 0o- On the other hand these same plots can be 
easily obtained from luminosity observed profiles. Therefore we can use the points in this 
plane to give a gross diagnostic of the models more appropriated to each class of objects. 
Observe that all models tend to converge to a single point for low —00 models, independent 
of the q index. This is a consequence of the density profile as a function of the potential as 
shown in figure [H For low values of — 0o the profiles tend to approximate a polytropic with 
n = 5/2. As a consequence, these models tend to present a unique surface density profile 
and therefore also have a unique pair of concentration indices Ci ~ 1.46 and C2 — 1.28. This 
regime coincides with the profiles of highly tidally affected models where the tidal radius 
is closer to the core radius. Therefore, strongly tidally truncated GC tend to be equally 
well represented by different q models. However, when we begin to sample larger values 
of the central potential we observe the differences among the q models. In this regime the 
individual density profiles behaves quite differently as a function of the parameter q. 

The King models, represented here by the subclass with q = 1, are known to give a good 
representation of the globular clusters. An empirical fit to the radial distribution of GC 



show that [391] 



pW=Po(1 + ^)-" (19) 

with 3.5 < a < 4. Using this empirical description we indicate in figure H] the region 
corresponding to the GC profiles by the letter G. We can see that models having a = 3.5 



tend to fall quite close to the region described by the King models. At this point one may 
ask if these q-lowered models may provide at least a reasonable description for elliptical 
galaxies. Such a possibility will be examined next section. 

IV. PROFILES FOR ELLIPTICAL GALAXIES 

The density profile of elliptical galaxies have been widely investigated either from the 
observational and theoretical point of view. A very successful approach consist in represent 
the luminosity profile of these objects by r^/"^, usually termed de Vaucouleurs law 
Although existing many indications that such empirical relation is not universal, there is a 
firm believe that it provides a good representation for several elliptical galaxies. 

A generalization of the de Vaucouleurs law is the so-called Sersic law where the luminosity 
profiles are represented by an expression, 

h = /elO-''"[(^)''"-^l (20) 

where n is the Sersic index and 6„ is a properly defined constant so that the effective radius 
Te contains half of the total luminosity. In particular, for n = 4 the de Vaucouleurs profile 
is recovered. Models with ^ 3 can be represented by the lowered King models, but those 
with n ~ 4 are definitely better represented by a model with 0.9 < g < 1. 

A simple integration of the de Vaucouleurs law leads to profile corresponding to a pair of 
concentrations indexes Ci = 2.730 and C2 = 2.506, which is represented by the label r^^'^ in 
figure m It is interesting to observe that this point is closely represented by a model q ~ 0.95 
and —00 — —10.45. in figure [5] we present the projected density profile corresponding to 
this model. As one may see, this model is a good representation of the de Vaucouleurs law 
over an interval of nearly 10 orders of magnitudes. There are two points that deserve a 
further comment. The first one is that models with parameters q =~ 0.95 and = —10.45 
are the best representation of the de Vaucouleurs law. If we explore the parameter region 
neighboring this point we can verify that the agreement is progressively lost. In figure H] we 
illustrate that point by exploring points departing 1% from that optimum model. 

The other point of interest is that our representation of the de Vaucouleurs law does 
not exactly reproduce the r^/^ profile. There is a noticeable oscillation around that trend. 
Over the 10 magnitude interval were the q profile is close to the de Vaucouleurs law this 
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FIG. 3: Projected density profiles for the 3 classes of models presented in Figure [2j Models 
corresponding to lower values of —(j)o, and hence C, have a higher similarity due to the dominant 
effect of the external polytropic structure (n = 5/2). 

oscillations have an amplitude of the order of 0.2 magnitudes. The remarkable fact is that 



the same phenomenon was observed in 



38l | in the elliptical galaxies belonging to the Virgo 



cluster. This is clearly a point that deserve a closer scrutiny in the future. 



V. CONCLUSION 



In this paper we have proposed a new family of lowered models with basis on the nonex- 
tensive kinetic theory. As we have seen, these q-lowered models based on Tsallis statistics 
extend naturally the class of Michie-King family and are also able to reproduce the observed 
structure of globular clusters. Although giving a broader range of density profiles than the 
usual King-Michie lowered models, the region occupied by these system in the concentration 
index diagram remains close to the q = 1 models which correspond exactly to the classical 
King models. It remains a matter of investigation to verify if the actual observed profiles of 
galactic globular clusters do favor a g = 1 lowered Maxwellian distribution or if there exist 
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FIG. 4: This concentration diagram shows the gross characteristic of the mass distribution as a 
response in changing the q nonextensive parameter. The interrupted Une with the label G mark 
the approximated location of GC lying close to the King model, corresponding to g = 1. The 
circular region with the label r^/"^ law marks the position of the elliptical galaxies obeying the de 
Vaucouleurs law. In this case, a good description of the luminosity profile is provided by a model 
with q = 0.95. 

room do discuss the presence of departures from this solution. 

For elliptical galaxies the situation is much more favorable to the nonextensive models. 
It is known that the Michie-King models do not give a close representation of the observed 
profiles of elliptical galaxies. Some objects, like NGC 4472, possibly affected by tidal trun- 



cation are well represented by the King profile {4]. However, the vast majority of bright 
elliptical, as NGC 3379, do no fit quite well in the King model and tend to be closer to the 
de Vaucouleuers profile [9j]. In this concern, the contribution of the lowered nonextensive 
models can be quite important since it is possible to find models which are closely resembling 
the structure of elliptical galaxies. In fact a model with q = 0.95 m is able to reproduce 
the de Vaucouleurs empirical law over more than 10 magnitudes of surface brightness. It 




FIG. 5: A comparison of the q = 0.95 model with the de Vaucouleurs law showing that this model 
may give a good description of the r^/^ law over an interval of 10 magnitudes. 

could be a matter of great interest to pursue a more detailed comparison of these models 
with the actual data on elliptical galaxies. In particular these non extensive models predict 
a systematic departure from the de Vaucouleurs law in the central region similar to the find- 
ings by Ferrarese et. al Using data from the Hubble space telescope these autors have 
demonstrated that the very internal regions of ellipticals tend to show an internal isother- 
mal core quite different from the prediction of a pure de Vaucouleurs law. Using bright 
ellipticals in the Virgo cluster they show that the brightness profile become closer to the de 
Vaucouleurs law only after the inner 10 arcsecond region. There is also a prediction of an 
external truncation but this could be a more difficult task to detect since this truncation is 
heavily dependent of the sky subtraction procedure. Finally, we would like to stress that the 
unified treatment for globular clusters and elliptical galaxies based on Tsallis distribution as 
presented here may, in principle, be extended to the framework of the so-called Kaniadakis 



non-gaussian k-distributions 
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